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ABSTRACT 



Sun-like stars show intensity fluctuations on a number of time scales due to various 
physical phenomena on their surfaces. These phenomena can convincingly be studied 
in the frequency spectra of these stars - while the strongest signatures usually originate 
from spots, granulation and p-mode oscillations, it has also been suggested that the 
frequency spectrum of the Sun contains a signature of taeniae. 

We have analyzed three stars observed for 13 months in short cadence (58.84 seconds 
sampling) by the Kepler mission. The frequency spectra of all three stars, as for the Sun, 
contain signatures that we can attribute to granulation, taeniae, and p-mode oscillations. 

The temporal variability of the signatures attributed to granulation, taeniae and 
p-mode oscillations were analyzed and the analysis indicates a periodic variability in 
the granulation and taeniae signatures ~ comparable to what is seen in the Sun. 

Subject headings: stars: activity — stars: individual(KIC 6603624, KIC 6933899, KIC 
11244118) — stars: oscillations — stars: solar-type 



1. Introduction 



Low-mass main sequence stars (Sun-like stars) have an outer convection zone, which means 
that they show granulation on their surface. The outer convection zone is responsible for exciting 
acoustic oscillations inside these stars, which, at least i n the stars with ma sses close to the Sun, 
have observable amplitudes at the surfaces of these stars (jHoudek et al.lll999l ) . When observing the 
disk-intergrated intensity of these stars as a function of time, the resulting light curves will thus 
include the signatures of both granulation and oscillations. These signatures are usually hard to 
identify directly in the light curves as granulation, which has the largest amplitude, is a strongly 
non-coherent phenomenon. The same applies to oscillations in Sun-like st ars as their oscillation s 
are damped, which means that they are also not coherent on long timescales (iGoldreich et al.lll994l ). 
An exception from this are evolv ed stars, which can have clear signatures of both granulation and 
oscillations in their light curves (ISchwarzschildlll975l ). 



Granulation and oscillations are therefore usually studied in Sun-like stars by analyzing the 
frequency spectrum of their light curves. These frequency spectra show a characteristic decline in 
power with increasing frequency caused by granulation. The location in frequency of this charac- 
teristic decline gives the characteristic time scale of granulation - equ ivalent to the life time or the 
turn-over time of t he granules, which has been measured in the Sun (lHarveylll985l ). sun-like stars 



(jMichel et al.ll2008l ) and evolved red-giant stars (IMathur et al 



2011 
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The characteristic frequency of the p-mode oscillations in Sun-like (and evolved) stars (i.e. 
the frequency of maximum power) is proportional to the atmospheric acoustic cut-off frequency 
which is given as the sound speed divi ded by twice the pres sure scale height near the surface - 



assuming the atmosphere is isothermal ([Belkacem et al.ll201ll ). Traditionally, it has been assumed 



that the chara cteristic time scale of granulation would also scale wi th the atmospheric acoustic cut 



off frequency ()Brown et al 



1991 



Kjeldsen fc Beddind 1199,1 



2011 



This means that there would 
be a linear relation between the atmospheric acoustic cut-off frequency and the difference between 
the characteristic time scale of granulation and the p-mode oscillations. 

It has been known for some time that the frequency spectrum of the Sun shows another feature 
in the region located between the charac teristic decline in p ower caused by granulation and the 
bump caused by the p-mode oscillations (iHarvey et al.lll993l ). Lately, it has been p roposed that 
this feature is likely caused by bright faculae on the surface of the Sun (lKaroffll2012l ). 



At l ow frequencies, pe aks can be seen caused by the rotational modulation of long-lived 
sunspots (jLanza et al.ll2004l ). At frequencies below 1 //Hz a linear decline in power with increasing 
frequency ca n be seen, which is caused by the weak coherency of activity-related phenomena such 
as sunspots (|Harvev et al.lll993l ). At frequencies higher than the atmospheric acoustic cut-off fre- 
quency evidence is seen for the so called high-frequency peaks, which could be traveling waves or 



chromospheric osc illations () Garcia et al 



2008 



Karoil2009|). 



1998 



Jimenez et al 



2005 



Karoff 



2007 



Karoff &: Kjeldsen 



At least two different interpretations of the physical processes responsible for the different 
features in the solar freque ncy spectrum can be f o und in the literatur e . In the first interp reta- 
tion, which has bee n used bv I Andersen et al.l (|l998l ). lAnklin et al.l (|l998l ). lAigrain et al.l (|2004l ) and 
Michel et al.l (|2009l ). the different features are explained by: rotation, supergranulation, mesogran- 
ulation, granulation and osc il lations (see Table 1). In the seco nd inte r pretat ion, which has been 
used by iHarvey et al.l (jl993l ) , I Vazquez Ramio et al.l (j2002l ) and iKarofB (j2012l ) , they are explained 
by: rotation, granulation, bright points (or faculae) and oscillations (see Table 2). Though this 
is a rough separation of all these studies into only two different interpretations, it exemplifies the 
general difference between these studies - that the studies using the first interpretation ascribe the 
feature just below 1000 /iHz to mesogranulation and the feature just above to granulation, whereas 
the studies using the second interpretation ascribe the feature just below 1000 /xHz to granulation 
and the feature just above to bright points or faculae. 



Following iRoudier et al.l (jl998l ) the three different convective scales can be identified as: Su- 
pergranulation characterized by a mean size of 20 to 50 Mm, mesogranulation by a mean size of 
5 to 10 Mm and granulation by a niean size of aro und 1000 km. The coherence time-scale of 
granulation is bet ween 2 and 6 min. (Del Mordl2004l. and refe r ences herei n). This agree s with the 
interpretation by iHarvev et al.l (|l993l ) , IVazquez Ramio et al.l (j2002l ) , and iKarofa (j2012l ) that the 
signature just below 1000 fiHz is due to granulation. Mesogranulation and supergranulation would 
have to have significant longer coherence time-scales and thus manifest themselves at frequencies 



- 4 - 



significantly below 10 00 //Hz. No such manifestations seem to be visible in the observed frequency 
spectrum of the Sun (|Karofj|2012l ). 



That the manifestation of mesogranulation and supergranulation is not observable in the fre- 
quency spectrum of the Sun does not mean that mesogranulation and supergranulation does not 
contribute to the background in the frequency spectrum. The contribution is just not large enough 
to be observable. The same is tr ue for the second population of granules wi th cohe r ence t ime-scales 
of around 1 min. identified bv iDel Mord (j2004l ) . It has been shown by |Karoflg (j2012l ) that this 
second population of granules, which would be both more numerous and with a lower brightness 
contrasts, would not be observable in the frequency spectrum, as the amplitude of the signature in 
the frequency spectrum scales with the brightness contrast over the square root of the number of 
granules. We thus follow the second interpretation in this paper. 

Here we analyze observations from the Kepler spacecraft of the three Sun-like stars: KIC 
6603624, KIC 6933899 and KIC 11244118. KIC 6603624 was o ne of the three stars presented in 
the first results on solar- like stars based on Kepler observations (jChaplin et al.ll2010l ) and already 
in that work it was noted that the frequency spectrum of this star shows a number of different 
components. One component possible originating from bright facul ae (as seen iii the S un). The 
eigenmode frequencies of all three stars were analyzed in detail by iMathur et al.l (l2012l ) and the 
resulting stellar parameters are shown in Table 3. Note that the uncertainties in Table 3 are formal 
uncertainties only, which do not include estimates of systematic contributions. 

The three star s were selected partly because they were among the 22 stars modelled by 
Mathur et al.l (j2012l ). which provides us with precise stellar parameters from asteroseismology and 
partly because they all show frequency spectra with features close to those observed in the Sun. 
The three stars are not the only Sun-like stars observed by Kepler that show the signature that we 
attribute to faculae in their spectra. 

The outline of this paper is the following: In section 2 we describe the 13 months of observations 
from the Kepler spacecraft that were used in this study and how these observations were pre- 
processed. A description of the analysis of the frequency spectra is given in section 3, including 
how we measure the parameters of the components that we attribute to granulation, faculae and 
p-mode oscillations. The results are presented in section 4 and conclusion are given in section 5. 



Observations 



The three stars were all observed in short cadence (58.85 sec. sampling. ICilliland et al.ll2010l ) by 



the Kepler mission (jKoch et al.ll2010). The specific obser ving times are give i i in T able 4. The light 



curves were prepared as described in I Jenkins et al.l (120101 ) and lGarcia et al.l ()201ll ). which removes 
most instrumental artefacts. The one artefact that was left in the light curves was a general declining 
trend. As all light curves are normalized to one, it is thus problematic to connect different quarters 
of observations. This is partly why we have chosen to calculate individual frequency spectra of 



- 5 - 



all months and then combining these spectra to an average spectrum ~ and because we do not 
need high-frequency resolution for studying granulation, faculae and the envelope of the p-mode 
oscillations. 

Apart from the declining trend, all three stars show RMS scatter with an amplitude of around 
one hundred ppm in the time domain and the noise level is relatively constant from month to 
month. The first month in quarter 2 contains a gap, as does the last month of quarter 7 and the 
second month of quarter 8. Nothing is done to these gaps in the analysis as the normalization to 
power density (see below) accounts for the missing observations. This means that the measured 
parameters are not affected by the gaps. 



3. Analysis 



The analysis of the signatures attributed to granulation, faculae, and p-mode oscillations in 
the frequency spectra consists of three steps. First, the frequency spectra are calculated from the 
light curves, then secondly, the resulting observed spectra are compared to a model in order to find 
the best model and thirdly, the parameters of the components attributed to granulation, faculae 
and p-mode oscillation are measured. 



The power density spectra were calculated using the least-squares method (jLomblll976l : iKarofi 



20081 ). Each spectrum was normalized by the effective observation length given as the reciprocal of 
the area under the window function in order to convert the spectra into power density and insure 
that the spectra obey Parseval's theorem. 

We have only analyzed the spe ctra in the f requency range between 100 /^Hz and up to the 
frequency of maximum power (as in lKaroffll2012l ). The reason for the lower limit is that this part 
of the spectrum is often affected by instrumental and rotational effect s, which are n ot the subject 
of this study. The reason for the higher limit is discussed in details by iKarofj ()2012l ) and is related 
to the choice of model, as we discuss below. 



3.1. The model 



Traditionally, the s ignatures attri buted to granulation and faculae have been evaluated using 
the model suggested bv iHarvevI ()l985l ): 



1 + {2ttvtY' 



(1) 



where f{v) is the power density at frequency v, a is the amplitude of the signature and r is the 
characteristic time scale. 



As discussed in iKarofj (j2012l ) various corrections have been made to this equation in order to 
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be able to model the frequency spectrum of the Sun. Here we adopt the model from iKarofla (|2012l ): 



Cgran'^gran'^gran 
1 + (27ri/rgran)^-^ 



+ 



CfacCfac7"fa 



+i?oscexp 



1 + (2^Z.Tfec)6-2 



+ Noise. 



(2) 



We have included the normalization constants Cgran and Cfaci which we will discuss in the next 
section, .ffosc is the height of the p-mode oscillation envelope, z^max is the frequency of maximum 
power and w is the width of the p-mode oscillation envelope. A white-noise term has also been 
added to the model. 

For the three stars in this study we adopt the s olar value of the exponents in eq. 2 - i.e. 
—3.5 lb 0.3 for granulation and —6.2 ± 0.7 for faculae (jKaroffll2012l ). The size of the exponents is 
a measurement of the amount of memory in the physical process responsible for the component. 
A larger exponent means less memory in the process. As the physical processes responsible for 
granulation and faculae are expected to be the same on other Sun-like stars, with global parameters 
close to the solar ones, the two exponents are expected to be close to solar values. 



As discussed by iKarofj (j2012l ) the model in eq. 1 fails to reproduce the observed solar acoustic 
background for frequencies higher than the atmospheric acoustic cut -off frequency. The re ason for 
this is that granulation ca nnot be n iodell ed with turbulent cascades (INordlund et al.lll997l ) as it is 
done in the drift model by iHarveyi (119851 ) . Turbulence shows a distribution with a slope of around 
— 2 in power, convection on the other hand has a lower limit in the time domain on which changes 
can take place. This means that on small time-scales (or at high frequency) convection is not 
noisy whereas turbulence is. Therefore in order to model the observed solar acoustic background 
for frequencies higher than the atmosp heric acoustic cut-off frequency an extra term has to be 
included in eq. 1 (or 2) as it was done in iKarofj (j2008l ). In order not to have to make assumptions 
about the behavior of this extra term in other stars, we have chosen not to model the high frequency 
part of the spectra and this is why we only analyze the spectra up to the frequency of maximum 
power. 



3.2. Normalization 

The amplitudes of the components attributed to granulation and faculae in eq. 2 can be 
understood as the variance these signals have in the time series. This is insured by Parseval's 
theorem, but in order for eq. 2 to obey Parseval's theorem, eq. 2 needs to be normalized using the 
normalization constants Cgran and Que- 

Parseval's theorem says that: 

'^^'^Ny t=tN I /,N -|2 

E /(-)^- = E ^^^v^^^ = -^ (3) 

i/=0 t=to 
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where v^^y is the Nyquist frequency, N the number of observations, o{t) is the observed intensity 
at time t and o the mean value of the observations. 5t is the time step and 5v is the frequency 
resolution {5v = j^)- The first sum should be summed over all natural frequencies, whereas the 
second sum should be summed over all time steps. As we assume the different contributions to the 
frequency spectra - granulation, faculae and the p-mode oscillation envelope, to be non-coherent 
and uncorrelated, Parseval's theorem can also be written as a sum of k different contributions: 



Ny 



5: Y.M^)su=a^=Y: 



-I 



(4) 



u=0 



and the different contributions can thus be normalized individually to Parseval's theorem, as they 
are additive in power: 



'-Ny 



(5) 



i/=0 



Following iMichel et al.l (|2009l ) it is thus possible to calculate the normalization constants Cgran and 
Cfac from eq. 2: 



'-Ny 



E 

u=0 



1 + (27ri/rJ"fe 



6u = al, 



(6) 



where are the exponents in eq. 2 ~ i.e. 3.47 for granulation and 6.20 for faculae. Assuming that 
i^j\[yT ^ 1 it can be shown that: 

(7) 



Using this formulation we obtain: Q 



Ck = 2Qfcsin 

gran = 5.46 and Cfac 



ak 



6.02. 



3.3. Minimization 



The model described above wa s matched to the observed frequency spectra using maximum- 
likelihood estimators as described in lKarofa (|2012l ) by calculating the logarithmic likelihood function 
i between N independent measurements Xi - i.e. power density at a given frequency, and the model 
fi given by a set of parameters A: 



N 



^In p{xi,X), 



(8) 



i=l 



where p{xi, A) is the probability density function, which is obtained from I Appourchauxl (j2004l ): 



p[Si{x,n),fi{X)] 



(9) 
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where Si{x,n) is the observed spectrum binned over n bins, F is the Gamma function and ^Uj and 
Vi are given as: 

k=i+n/2 

{n+l) E /'^W 

k=i—n/2 



and 



k=i+n/2 
fc=j— n/2 
'A;=i+n/2 

k=i—n/2 



(10) 



(11) 



fc=i+n/2 

E /^'(^) 

A:=i— n/2 

Note that is different from the u used for the frequencies in the frequency spectra. 

The general idea is here to compare a binned version of the observed spectra Si{x,n) to 
the model in order to calculate the likelihood - instead of assuming that the differences between 
the observed frequency spectra and the model are given either by a normal distribution or a 
distribution with 2 degrees of freedom. The spectrum should be binned over so many bins that 
the model becomes a good representation of the observed spectrum. On the other hand n should 
not be so large that the different features in the frequency spectra cannot be observed. We have 
thus used n = 100. If n is close to one eq. 9 reduces to the well known case from the modeling of 
individual p-mode oscillation modes where the differences between the o bserved freque ncy spectra 
and the model are given by a distribution with 2 degrees of freedom ( Gabriellll994l ). 



The minimization o f the likelihood fu nction was done using the simplex method as applied in 
the AMOEBA function (|Press et al. l ll992l l. 



3.4. Uncertainties 



Formal uncertainties can be calculated as the diagonal elements of the inverse of the Hes- 
sian matrix, but such uncertainties are only internal uncertainties that says something about the 
curvature of the minimized likelihood function. 

The uncertainties also have to take into account the averaging that is performed in the mod- 
eling. First we average individual spectra to an average spectrum and then we bin this (and the 
individual spectra) over n bins. While the first effect can be acco unted for by dividin g the un- 
certainties by the square root of the number of individual spectra (jAppourchauxl l2003l ) it is not 
straight forward how to account for the second effect. 



If the differences between an observed frequency spectrum sampled at the natural frequencies 
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and the model were assumed to be given by a distribution with 2 degrees of freedom and this 
spectrum was binned over n bins, then we know that the differences between the binned spectrum 
and the model would be given by a distribution with 2n degrees of freedo m. In that case the 



uncertainties would have to be adjusted accordingly (jAppourchaux et al.ll2012l ). As the probability 
density function that we have used (eq. 9) accounts for the binning, it is, on the other hand, not 
clear if any corrections are need to the uncertainties. 

To avoid inconsistencies in the quoted uncertainties, we have chosen to give the average values 
from the results of different monthly, individual spectra rather than results of an average spectrum. 
In this way we can also give the uncertainties as the uncertainties on the mean value. This is on 
the other hand not possible to do on the temporal results, so here the error bars only represent 
formal uncertainties. 



4. Results 

Using the method described above we have modelled a summed spectrum calculated as the sum 
of individual spectra for each of the 13 months of observations for all three stars. These summed 
spectra were used to calculate the significance of the component attributed to faculae using the 
formulation described above. The average spectra are shown in Fig. 1 together with the model with 
(red line) and without (blue line) the component attributed to faculae. We have also analysed the 
individual spectra for each of the 13 months of observations for all three stars and the measured 
mean values and their temporal variability are shown in Figs. 4-7. 



4.1. Significance of the component attributed to faculae 

The significance of the signature of the component attributed to faculae was calculated using 
the logarithmic likelihood ratio A: 

In A = £(Ap+g)-£(Ap), (12) 

where p is the number of free parameters in the model without the component attributed to faculae 
(i.e. 6) and q is the number of additional parameters in the model with the component attributed 
to fa culae (i . e. 2), and comparing the value of — 21nA to a x^ distribution with q degrees of freedom 



as m 



Karoffl (12013). 



A visual inspection of the frequency spectra reveals that the model with the component at- 
tributed to faculae matches the observed spectra better than the model without such a component 
in all three cases. This is also reflected in the returned logarithmic likelihood ratios, which are 
shown in Table 5. The low ratios mean that the significance of the component attributed to facular 
is one in all three stars - i.e. the probability of having a logarithmic likelihood ratio smaller than 
e.g. -12 returned from a x^ distribution with 2 degrees of freedom is larger than 0.99999. 
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in order to test that the measured significances of the faculae components were not affected 
by using solar values for the exponents in eq. 2 we also calculated the logarithmic likelihood ratios 
allowing the exponents to vary around the solar values. This did change the returned logarithmic 
likelihood ratios, which are also shown in Table 5, but not the conclusion that the significance of 
the component attributed to faculae is one in all three stars. In other words it is clear that the 
conclusion that the component attributed to faculae is significantly present in all three stars is valid 
whether the exponents are fixed to solar values or allowed to change freely. 

When the exponents were fixed to solar values, the exponent of the granulation component 
were the same in both models, but when the exponents were allowed to change freely this was not 
the case - here the values of the granulation exponent were generally lower in the model with no 
facular component. Generally, for all the tests we did, the exponents were returned with values 
between 3 and 5 for the component attributed to granulation and 6 and 8 for the component 
attributed to faculae. For the model without the faculae component the exponents were returned 
with values between 1.7 and 4. 

The returned logarithmic likelihood ratios also agrees with the impression from a visual inspec- 
tion - that the signature of the component attributed to faculae is strongest in KIC 6603624 and 
weakest in KIC 6933899. In fact, as can be seen in Fig. 1, it was not possible to get a satisfactory 
agreement between the model with no component attributed to faculae and the observed spectrum 
of KIC 6603624. 



4.2. Mean values of the measured parameters 

The mean values of the measured parameters and their uncertainties were calculated from the 
values measured in the 13 spectra of the individual months (see Table 6). 

It is well known that there exist a relation between the frequency of maxii num power and th e 



large frequency separation, which again depends on the mean stellar density (jStello et al.ll2009l ). 
This agrees nicely with the fact that KIC 6603624 with a density almost twice as large as the two 
other stars is also the stars with the largest frequency of maximum power (see Table 6). The mean 
density of KIC 6603624 is still lower than the mean density of the Sun and so is the frequency of 
maximum power. 

Even though three stars is far from enough for calculating scaling relations we have tried 
to compare both the time scales and amplitudes of the various components to the frequency of 
maximum power in Figs. 2 & 3. It is seen that the time scales of both the components scale 
inversely with the frequency of maximum power (Fig. 2). We do not see any clear relation between 
the amplitudes of the componen t s and t he frequency of maximum power (Fig. 3). This contradicts 



the predictions bv lChaplin et al.l (|2011bl ): that the amplitudes should scale inversely with the square 
of the frequency of maximum power. The reason for this is most likely that we only analyze three 
stars in this study. 
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4.3. Temporal variability 

We have measured the temporal variability of the background and oscillation envelope param- 
eters in the 13 spectra of the individual months for each star. The temporal variability of the 
parameters of the components attributed to granulation and faculae are shown in Figs. 4-7. The 
white noise component is shown for comparison in Fig. 8. 

In order to test if any of the parameters show periodic variability we found the highest peak 
in the periodogram of the parameters as a function of time and performed a simple test of the 
significance of this peak. The test was performed by measuring the amplitude of the highest peak 
in one million artificially generated periodograms. These periodograms were generated by taking 
the uncertainties ascribed to each measurements at each time step (calculated as the diagonal 
elements of the inverse of the Hessian matrix) and multiply them with a normally-distributed 
random number with a mean of zero and a standard deviation of one. The significance was then 
calculated as the fraction of cases where the highest peak in these one million artificial generated 
periodograms had an amplitude lower than the amplitude of the highest peak in the observed 
periodogram. 

Here we adopt the definition that a period with a significance level above 0.99 is significant, 
whereas a period with a significance level above 0.95 is marginal significant. This leads to the 
conclusion that we see a significant period of 322 it 11 days (error bars calculated assuming normal 
distributed noise on the measured amplitudes) in the amplitude of the component attributed to 
granulation in KIC 6603624. Marginal significant periods are found in the amplitude of the compo- 
nent attributed to faculae in KIC 6603624 (257 it 11 days) and in the time scale of the component 
attributed to granulation in KIC 6933899 (233 it 13 days). The significant levels of the highest 
peaks in the periodograms of the temporal variability of the granulation and faculae parameters 
for the three stars are given in Figs. 4-7. No significant or marginal significant periods were found 
in the parameters of the envelope of the p-mode oscillations. 

The periods of the periodic variability seen in the components attributed to granulation and 
faculae in KIC 6603624 and KIC 6933899 are between two hundred and three hundred days. This 
could be equivalent to the quasi-annual modulation of the amplitude of the component attributed 
to granulation in the frequency spectrum of intensity observations of the Sun with the Variability 



of solar IRradiance and Gravity Oscilla tions (VIRG O jFrohlich et al.lll995l ) instrument on the Solar 



and Heliospheric Observatory (SOHO. lKaroflll20 1 2| ) . Unfortunately, it has been impossible to study 



such a modulation in veloc ity observations from the with Global Oscillations at Low Frequency 



(GOLF, iGabriel et al.lll995l ) instrument also on SOHO due to a variat ion of the observation height 



in the solar atmosphere induced by the orbital period of the satellite (jLefebvre et al.ll2008l ). 



Even though non of the identified periods are identical we did calculate the linear Pearson 
correlation coefficients between the granulation and faculae parameters (see Table 7) . The numerical 
values of these correlation coefficients were generally very low. The largest one being between the 
amplitude of the component attributed to faculae and the time scale of the component attributed 
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to granulation in KIC 6922899 with a value of 0.74. The low levels of these correlation coefficients 
verify that there are no correlations between the different parameters and thus that the identified 
periods are intrinsic to the amplitude of the component attributed to granulation in KIC 6603624, 
the amplitude of the component attributed to faculae in KIC 6603624 and the time scale attributed 
to granulation in KIC 6933899, respectively. 



4.4. Comparison of different analysis methods 



It was shown by iMathur et al.l (j201ll ) that different analysis methods with different free pa- 
rameters, models, data and number of components can provide different values of the estimated 
parameters, but that, for a given method, the results and the trends are consistent. Though we do 
not want to redo this analysis here it is still interesting to investigate if the component attributed 
to faculae can also be found using different analysis methods. 

We therefore redid the analysis of the frequ ency spectra of the three stars usirig the CAN 
and A2Z methods from iMathur et all koill ) (see iMathur et all boiol : iKallinger et^ boiol . for a 
detailed description of the methods) and a method that used the same model as described here, but 
no binning of the spectra and thus assuming that the differences between the observed frequency 
spectra and the model were given by a distribution with 2 degrees of freedom. 

The result of the comparison was that the component attributed to faculae was significantly 
present in the observed frequency spectra with a significance close to 1 in all three stars in all the 
tests where convergence to physically meaningful parameters could be reached. 

The comparison also showed that the uncertainties on the mean parameter values are realistic 
— i.e. they are comparable in size to the difference between the different analysis methods and 
the uncertainties returned by the MultiNest algorithm used in the CAN methods (which have been 
shown in brackets for comparison in Table 6). The error bars on the temporal results calculated 
from the diagonal element s of the inverse of the Hessian matrix are on the other hand generally 
too low (as also noted by Mathur et ablboill ). 



5. Conclusions 

The analysis of the observed frequency spectra of the three Sun-like stars KIC 6603624, KIC 
6933899 and KIC 11244118 from 13 months of high-precision, high-cadence photometric observa- 
tions from the Kepler spacecraft has revealed signatures of what is likely granulation, faculae, and 
p-mode oscillations in all three stars. 

It is seen that the characteristic frequencies (or equivalent, time scales) of the components 
attributed to both granulation and faculae inversely scale with the frequency of maximum power 
and that the frequency of maximum power scale with the mean density of the stars as expected as 
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the a tmospheric acoustic cut-ofF frequency scales with the large frequency separation (|Stello et al, 
20091 ). 



The analysis of the temporal variability of the measured parameters revealed periodic variabil- 
ity in the amplitude of the component attributed to granulation in KIC 6603624 with a significance 
level above 0.99 and a period of 322 ±11 days. Marginal significant periodic variability with a 
significance level above 0.95 were found in the amplitude of the component attributed to faculae 
in KIC 6603624 with a period of 257 it 11 days and in the time scale of the component attributed 
to granulation in KIC 6933899 with a period of 233 ± 13 days. 

The temporal variability in KIC 6603624 and KIC 6933899 could have a origin similar to the 
periodicities of around 1-2 yrs. seen in a nurn ber of indices related to solar act ivity, including: Sun- 



netic indices (Mursula et al 



spot number and neutri no flux (ISakurailll979l l. Galactic cosmi c-ray intensities (lyaldes-Galicia et al 



1996), flare o ccurrence ( Ant aloval 1994 ). solar wind velocities ( Richardsoii et al. 19941). 



20031 ) ■ coronal hole area and radio emission (Va. 



20081 ). the rotation of the Sun near the b ase of its convective zone (IHowe et al 



oo geomag- 



des-Galicia &: Velasco 



2000l) and lately in the 



residuals of the p-mode frequency shifts (jFletcher et al.ll2010l : iBroomhall et al. ISimoniello et al 
2OI2I ). 



The periods of the periodic variability in KIC 6603624 (and KIC 6933899) are between two 
hundred and three hundred days, which is somewhat lower than the 1-2 yrs. periodicities observed 
in the indices related to solar activity. This agrees with the lower densities of these stars compared 
to the Sun, but it is not clear if it agrees with the period of any possible dynamo in these stars. 
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Table 1: Interpretation of the physical processes resp o nsible for the differ e nt features in the so - 
lar frequency spect rum used bv lAndersen et alj ( 19981 ): lAnklin et al.l (|l998l ): lAigrain et alj ( 20041 ) : 



Michel et al.l (120091 ) 



0-1 //Hz 
10 - 100 fiRz 
80 - 1000 fiBz 
800 - 3000 nBz 
2000 - 4000 /iHz 



Rotation 
Supergranulation 
Mesogranulation 
Granulation 
Oscillations 



Table 2: Interpretation of t he physical processes responsible for the diff e rent f e atures i n the solar 
frequency spectrum used by iHarvey et alJ (jl993l ) ; IVazquez Ramio et al.l (|2002l ) ; iKarofla (|2012l ) 



- 1 fiRz 
100 - 1000 ^fHz 
1000 - 3000 /iHz 
2000 - 4000 //Hz 



Rotation 
Granulation 

Faculae 
Oscillations 



Table 3: Stellar parameters from asteroseismology (from lMathur et alJl2012l : IVerner et al.ll2011 



KIC 


6603624 


6933899 


11244118 


R{Rq) 


1.15 ± 0.01 


1.58 ± 0.01 


1.55 ± 0.01 


M(Mq) 


1.01 ± 0.01 


1.10 ± 0.01 


1.01 ± 0.01 


T,s (K) 


5416 


5616 


5507 


Kp (Mag) 


9.1 


9.6 


9.7 


Age{Gyr) 


8.51 ± 0.23 


6.28 ± 0.15 


8.93 ± 0.04 



Table 4: Observing times, 'a' means: KIC 6603624, 'b': KIC 6933899 and 'c': KIC 11244118. 



Quarter 


Start Data 








1 


13 May - 17 Jun 2009 


a 






2.1 


20 Jun - 19 Jul 2009 




b 


c 


5 


20 Mar - 24 Jun 2010 


a 


b 


c 


6 


24 Jun - 23 Sep 2010 


a 


b 


c 


7 


23 Sep 2010 - 6 Jan 2011 


a 


b 


c 


8 


6 Jan 2011 - 20 Mar 2011 


a 


b 


c 



Table 5: Logarithmic likelihood ratios [ln(A)] - with fixed and free exponents. 



KIC 


6603624 


6933899 


11244118 


ln(A)fixcd 


-531.67 


-420.3 


-98.35 


ln(A)fi.ee 


-683.16 


-21.34 


-53.34 
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Fig. 1.— Binned average spectra of KIC 6603624, KIC 6933899 and KIC 11244118 (from top to 
bottom). The solid red lines are the model with the component attributed to faculae and solid 
blue lines are the model without the component attributed to faculae. The dashed lines show the 
component attributed to granulation, whereas the dotted lines show the component attributed to 
faculae and the dash-dotted lines show the p-mode oscillations envelope. 
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Table 6: Mean parameters of the three stars and the Sun. Definitions of the different parameters are 
given by eq. 2. The number in brackets are the uncertainties returned by the MultiNest algorithm. 





KIC 6603624 




KIC 6933899 




KIC 1124418 




Sun 




^gran 


62.8 ± 1.5 


(1) 


74.5 ± 2.3 


(1) 


87.4 ± 2.9 


(1) 


62.4 ± 0.6 


ppm 


'^gran 


280.8 ± 5.9 


(4) 


474.3 ± 15.4 


(10) 


448.0 ± 10.7 


(3) 


214.3 ± 2.9 


sec 


CTfac 


76.5 ± 2.4 


(1) 


28.7 ± 7.5 


(1) 


74.7 ± 3.1 


(1) 


50.1 ± 0.1 


ppm 


Tfac 


66.1 ± 1.3 


(3) 


180.1 ± 2.4 


(7) 


153.1 ± 1.6 


(7) 


65.8 ± 0.3 


sec 


Hose 


1.4 ± 0.1 


(0.03) 


2.1 ± 0.2 


(0.09) 


4.5 ± 0.3 


(0.1) 


6.2 ± 2.5 


ppin^//iHz 


^jnax 


2477 ± 12 


(8) 


1368 ± 10 


(5) 


1426 ± 14 


(5) 


3104 ± 36 


/iHz 


W 


246.5 ± 0.9 


(10) 


138.9 ± 2.0 


(8) 


233.5 ± 2.9 


(7) 


316 ± 36 


fiHz 


Noise 


0.4 ± 0.1 


(0.1) 


4.1 ± 0.1 


(0.1) 


2.8 ± 0.4 


(0.1) 




ppm^//iHz 



^ 2 




0.4 0.5 0.6 0.7 0.8 0.9 
Frequency of maximum power relative to the Sun 



Fig. 2. — Measured time scales of the components attributed to granulation (red) and faculae (blue) 
as a function of the frequency of maximum power relative to the Sun. It is generally seen that both 
time scales decreases with increasing frequency of maximum power. 




0.4 0.5 0.6 0.7 0.8 0.9 
Frequency of maximum power relative to the Sun 



Fig. 3. — Measured amplitudes of the components attributed to granulation (red) and faculae 
(blue) as a function of the frequency of maximum power relative to the Sun. No general trend can 
be seen. 
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Fig. 4. — Temporal variability of the amplitude of the component attributed to granulation. Pe- 
riodic variability with a significance level above 0.99 and a period of 322 it 11 days is seen in KIC 
6603624. 
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Period = 257 days 
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Fig. 5. — Temporal variability of the amplitude of the component attributed to faculae. Marginal 
significant periodic variability with a significance level above 0.95 and a period of 257 it 11 days is 
seen in KIC 6603624. 
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Fig. 6. — Temporal variability of the granulation time scale. Marginal significant periodic variability 
with a significance level above 0.95 and a period of 233 ± 13 days is seen in KIC 6933899. 
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Fig. 7. — Temporal variability of the facular time scale. 

Table 7: Correlation between the faculae and granulation parameters (KIC 6603624 / KIC 6933899 
/ KIC 11244118) 





Cfac 


'^gran 


Tfac 




0.56 / 0.21 / 0.09 


0.15 / 0.26 / -0.47 


-0.72 / -0.56 / -0.61 


Cfac 




0.59 / 0.74 / 0.36 


-0.09 / 0.01 / -0.08 








0.06 / 0.18 / 0.18 



- 24 - 



S 2 





- KIC 6603624 




Period = 1 79 


days 


- p = 0.2185 








+ 
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Fig. 8. — Temporal variability of the white noise level. 



